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SUMMARY 

The motion of charged, conducting droplets present in an insulating fluid 
medium is analyzed under the action of an electric field, in microgravity. 
Previous analyses of this problem have considered the Maxwell stresses as the 
only driving force. In the present study, arguments from macroscopic thermody- 
namics and the molecular theory of surface tension are used to show that sur- 
face tension gradients can be induced due to the variation of the electric 
potential on the interface. In the limit of Reynolds numbers small compared to 
unity, the terminal velocity of migration of the droplet is calculated under 
the combined action of the Maxwell stresses and the surface tension gradients. 
The results show that there are no surface tension gradients (i.e., no electric 
potential variation at the interface) in a case that is considered strictly 
electrostatic. With the inclusion of surface currents due to the convection of 
the surface charges, surface tension gradients do exist and tend to reduce the 
terminal velocity of the droplet. The shape of the droplet altered by the 
motion has also been calculated, when the deformations from the spherical shape 
are smal 1 . 


NOMENCLATURE 

A^ area of the interface 

E electric field 

E 0 electric field far away from the droplet 

F e ,Fh net electric and hydrodynamic force on the droplet 

f deformed droplet shape 

J 1 , J2 constants defined by equations (53) and (54) 
kg Boltzmann's constant 

k£ electrical conductivity of the droplet 

N normal stress 

P pressure 

Q 0 ,q 0 dimensional total surface charge and dimensional average surface charge 
density of the droplet 
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q, q dimensional and dimensionless surface charge density 

R,r dimensional and dimensionless radial coordinate 
R a ,Rb principal radii of curvature 

Re Reynolds number 

R 0 droplet radius 

r 1 2 vector joining the center of mass of two molecules 

S entropy 

T temperature 

t unit vector tangent to droplet surface 

U, u dimensional and dimensionless radial velocity 

internal energy 

V , v dimensional and dimensionless tangential velocity 

Vr reference velocity 

y volume 

Voo.Voo dimensional and dimensionless terminal velocity of the droplet 

Y axial coordinate 

(3 tagged constant (equals 1 or 2) 

r, T 2 one and two particle (molecules) surface excess quantities 

Y tagged constant (equals 0 or 1 ) 

c electrical permittivity 

n dimensionless droplet deformation 

9 tangential coordinate 

X ratio of lengths of axes of the droplet in the flow direction to that 
in the normal direction 

y dynamic viscosity 

v kinematic viscosity 

a surface tension 
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-t,x e hydrodynamic and Maxwell shear stress 

4> azimuthal coordinate 

$,<j> dimensional and dimensionless electric potential 

particle (molecular) potential energy from many particle Hamiltonian 
dimensional and dimensionless stream function 
Subscripts : 

1 outer fluid 

2 droplet medium 

i interface 

o reference values 


INTRODUCTION 

Electrically charged droplets present in a fluid with which they are 
immiscible, will move when they are subjected to an electric field. They 
attain a steady speed of electrohydrodynamic migration, termed the terminal 
velocity, in a uniform field, when the electrical force is balanced by the vis- 
cous resistance. Such migration is also called electrokinetic or electropho- 
retic motion in colloidal systems, where the motion is due to the interaction 
of the applied filed with the electric double layer present at the interface 
and has been extensively investigated (refs. 1 to 6). Electrohydrodynamic 
migration, like thermocapillary migration, is useful in many applications in 
microgravity, where, in the absence of buoyancy, the locations of droplets and 
bubbles need to be controlled by various mechanisms. 

There have been few investigations of the electrokinetics of fluid-fluid 
interfaces (ref. 6). Most studies (refs. 2 to 6) consider the Maxwell stresses 
as the driving force for the flow. Levich (ref. 1) has analyzed the motion of 
mercury drops in electrolytes, representing the driving force as a surface ten- 
sion gradient, arising from the electrocapillary effect. It is the contention 
of the authors that the Maxwell stresses and surface tension gradients due to 
electrocapillarity, though each occurs because of the applied field, are inde- 
pendent driving mechanisms as far as the flow is concerned and hence can act 
simultaneously. The combined effect of the two is analyzed in this study for 
a model problem of the migration of a charged, conducting droplet in an elec- 
trically insulating medium. Since the droplet medium is conducting, the charge 
resides as a surface charge at the interface. This situation is the same as 
the "charge monolayer" problem considered by Spertell and Saville (ref. 4). 

In what follows, the thermodynamics of the system are considered to show 
that surface tension gradients can be generated by the variation of the elec- 
tric potential at the interface, and reasoning from a molecular basis is used 
to explain the origin of the surface tension gradient. The electrohydrodynamic 
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formulation for the migration of a spherical droplet is then presented and the 
potential and flow fields and the terminal velocity are calculated. The calcu- 
lated fields are lastly used to determine the slightly perturbed shape of the 
droplet. 


THERMODYNAMICS AT THE INTERFACE 

Presented below is a simplified view of the thermodynamics of the inter- 
face where the occurrence of a charge double layer is excluded. Consider a 
system consisting of an electrically insulating fluid with an electrically con- 
ducting immiscible droplet inside, under microgravity (fig. 1). Let the system 
change from one equilibrium state to another, during which the system can 
exchange heat and electrical charge with the surroundings, but no mass enters 
or leaves the system. Denoting the insulating fluid, the droplet medium and 
the interface region by the subscripts 1, 2 and i respectively and assuming 
that each subsystem (l,2,i) remains in equilibrium at the initial and final 
states for a reversible process with the droplet and the interface not trans- 
ferring charge to the insulating fluid, the change in the internal energy of 
the system is 


d & s = d + 

where 

d&i « T i dS i 

d^ 2 = T 2 ds 2 

d U\ * TidSi + 


111 + d 

(i) 

- P]d r } 

(2) 

- p 2 d 

(3) 

dA ^ + <}> j dQi 

(4) 


a is the surface tension (or interfacial tension; the two words are used 
interchangeably), Qj is the surface charge brought to the interface from 
infinity where the potential is zero and <j>. is the electric potential at the 
interface. Keeping the Intensive variables fixed, if the extensive variables 
are permitted to undergo a finite change, (ref. 7) the following expression 
may be written for U\ 


^ \ — T^ Si + orAi + 4>.Qi (5) 

Differentiating equation (5) and subtracting equation (4) yields the Gibbs- 
Duhem equation for the interface as 

Aida + Qi d$ i + S i dT i * 0 (6) 


If the surface is isothermal during the process, dT i = 0 and hence 

°i - 

da = - T 1 d<j> - (7) 

A i 1 

which is analogous to the Lippman equation for electrocapillary phenomena 
(ref. 8). From equation (7), we may write 

V s a = -qi V s <j>. (8) 
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where is the surface gradient operator and q. is the surface charge den- 

sity. Thus surface tension gradients can occur when the surface is charged 
and the potential varies on the interface. 


MOLECULAR BASIS OF SURFACE TENSION GRADIENT 

If one thinks of the surface of the charged droplet as a membrane with 
embedded charge, then whenever there is a nonuniform charge density, there will 
be unbalanced electric forces between the charges. This is the physical pic- 
ture that persuades us that if the externally applied electric field induces a 
nonuniform charge, then a surface stress is induced. 

Referring to a treatise on the molecular theory of surface tension, one 
can write a molecular level expression for a as (ref. 9) 

a = -Tk-gT + ^ 

where T_ and T 2 are the 1 and 2 particle surface excess quantities respec- 
tively; 4> is the particle potential energy from the many particle Hamiltonian; 
i is a unit vector tangent to the surface; and kgT is the Boltzmann constant 
times temperature. The potential energy term above effectively accumulates the 
surface component of the force acting on a particle by all the other particles. 
T 2 ensures that contributions occur only where there is an excess or depletion 
of particle pairs with respect to the bulk. 

On our charged surface, $ will contain contributions from the external 
electric field, which is a solution to the electrostatic problem, as well as 
the local electric field from neighboring charged particles. It is the local 
electric field present in the nonuniform surface charge distribution case that 
we contend physically mandates a surface tension gradient in addition to Max- 
well stresses calculated from an external electric field. 


f 2 (? ? 


12 


$) • t 


dr 


12 


(9) 


ELECTROHYDRODYNAMIC FORMULATION 

Consider a conducting droplet present in an insulating fluid of infinite 
extent (Fig. 1). The droplet carries a total surface charge, Qp, and is sub- 
jected to a uniform electric field E 0 . We would like to calculate the termi- 
nal velocity of the droplet and the flow and potential fields in both fluids. 
As the droplet moves at its constant terminal velocity, it is convenient to 
choose a coordinate system on the droplet with the origin at its center of 
mass. In this coordinate system, the outer fluid approaches the droplet with 
the terminal velocity V , which is to be determined. 

The flow is considered to be incompressible, laminar and with spatially 
constant physical properties (viscosity, electrical conductivity, etc). The 
droplet is assumed to be spherical. From the geometry and the boundary condi- 
tions, the problem is symmetric about the flow direction (the droplet polar 
axis). The basic equations for the flow in the two fluids are the Navier- 
Stokes equations. The treatment of electrical effects in this study is the 
same as in the study by Melcher and Taylor (ref. 10). As mentioned in 
reference 10, the dynamic currents are so small and the induced magnetic field 


is negligible and hence, the basic equations are the electrostatic form of Max- 
well's equations. It is assumed that the current density is proportional to 
the electric field in the droplet (Ohms law), with the electrical conductivity 
being a constant. In both fluids, the electric displacement is taken to be 
proportional to the electric field, with the permittivity being a constant. 

The Maxwell equations then reduce to an equation for the electrostatic poten- 
tial. Subscripts 1 and 2 are used to denote the outer fluid and the droplet 
medium respectively. The coordinate system is R ,0 ,<t> with the origin at the 
center of mass of the droplet (fig. 1). 9 is measured clockwise from the 

point of incidence of the flow. Rationalized MKS units are used for the elec- 
trical quantities. 

A reference velocity scale Vr for the flow field may be estimated from 
the fact that when the droplet moves at the terminal velocity, the Coulombic 
force on it is balanced by the viscous forces. Hence, 

**0 Vo ~ 4 * R 0 2< “1 V R /R o ,=>V R - < l0) 

where q 0 = Q 0 /(4irR 0 2 ) is the surface charge density. Vr thus yields typical 
values of the velocities thay may be anticipated. The Reynolds number is 

Re = V D R /v = q^E R 2 /(}i 1 v) (11) 

R 0 0 0 0 1 


Re compares the inertial forces on a fluid element to the viscous forces on 
it. In what follows, Re for both fluids is assumed to be small compared to 
unity. The Navier-Stokes equations are written In the Re << 1 limit for the 
stream function (ref. 11). Dimensionless variables are defined as follows 


r 




v 





E o R o 


(12) 


where Vr is as defined by equation (10) and the stream function ^ is such 
that 


i at 

r 2 sln 6 88 


V = 


1 8^ 
r sin 0 9r 


The basic equations in the two fluids are 

d\ = 0 


where 


V 4> = 0 


D = 


9r‘ 


i_ e . n A a_ r i a_1 

2 s 80 Lsin 0 80 J 


v 2 - 1- i- 

v 2 9r 
r 


r 2 3_' 
8r 


r sin 0 


fe [ sin 9 y 


(13) 

(14) 


(15) 
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The boundary conditions are r -*■ ® are 


00 2 ■ 2 n 
-> 2 — r sin 0 


(16) 


i.e., the velocity tends to the free stream velocity v^ = VJV r and 

4> -| -r cos 0 (17) 

i.e., the electric field is uniform far away from the droplet (E* = -V$ = -Eq^). 
The boundary conditions at r = 1 (the droplet surface) are 
(a) = constant = 0, = constant = 0 (18) 

i.e., the radial velocities are zero at the surface. The constants may be set 
equal to zero without loss of generality. 


(b) 


3t 1 3l»2 

W~ = 3r~ 


(19) 


i.e., the tangential velocities are continuous, 


(c) 


3<t>l 3<l>2 

30~ ~ 30~ 


(20) 


i.e., the tangential fields are continuous, 


(d) 


q = 


3$i 


9$, 


-e. 


3R + e 2 3R J 


( 21 ) 


R=R. 


This boundary condition is written dimensionally and represents the discontinu- 
ity in the normal displacement fields in proportion to the surface charge den- 
sity (ref. 12). Since the total charge on the droplet is Q 0 = 4irR 0 2q 0 , we 
must have 


t r 


q 2irR^ sin 0 d0 = 4irR^ 2 q^ 
H 0 0 0 


( 22 ) 


(e) Shear stress balance: The balance of shear stresses at the interface 

may be written as 


e _ e ]_ 9o_ _ ^ 
t 1R© " T 2R0 + t 1R© " t 2R0 + R q 30 


(23) 


where = pR 


3R 


R=R 


is the hydrodynamic shear stress and 


t - cE E 
R0 R 0 


is the electric or Maxwell shear stress. From equation (8), 


R=R 
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3o/ 30 = - q 3<j>/30. Using dimensionless variables and defining q = q/q 0 , the 
shear stress balance at r = 1 may be written as 


3 

f V ll 

^2 

3 

y _2 

3r 

_r 

Pi 

3r 

.r . 



(24) 


It may be noted that - t 2 ®q = 3a/30 = -q< 3<(>/ 39) at r = 1 and hence the 

coefficient in the right hand side of equation (24) is 2. Since this is the 
only place in the formulation that the surface tension gradients are encoun- 
tered (surface tension is also encountered in the normal stress balance at the 
droplet surface and will be addressed later), the term 2q3<t>/39 will be tagged 
for purposes of comparison and written as (Jq3<j>/39. Hence (3 = 2 wi 1 1 recover 
the present formulation while |3 = 1 represents a formulation where either Max- 
well shear stresses or surface tension gradients (but not both) are considered. 


(f) Charge conservation: Since the interface can move with respect to the 

bulk of the droplet, the motion of the surface charges constitutes an electric 
current at the surface. Since charge must be conserved, and the outside medium 
is insulating, there must exist a radial conduction current at the surface from 
the inside medium. Using Ohm's law (current density = conductivity x field) 
and inventorying charge movement, charge conservation may be written as 


34> 2 

w 



sin 0 


I (<,v 


sin 9) 


(25) 


(g) Net force on the droplet: As the droplet moves at its terminal veloc- 

ity, the net force acting on it must be zero. Since surface tension is an 
internal force for the droplet and cannot act on itself, the sum of the hydro- 
dynamic and elertricaj forces exerted on the droplet by the outside medium must 
be zero, i.e., T e + Tf, = 0. By symmetry, the net force perpendicular to the 
flow direction is zero. That the net force parallel to the flow direction is 
zero may be expressed as 


irlr 


L C l E Rl E 01 S ’ n 6 2 E R1 C0S e J + 


h R ok 

X (jAin e - (zm, p 1 - p, 


cos 9 


sin 9 d0 = 0 


R=R 


(26) 


where Pi is the hydrodynamic pressure in the outer fluid. This condition is 
essentially what determines the unknown eigenvalue V^. Instead of integrating 

the surface stresses, F e may also be calculated alternatively by recognizing 
that it is equal to the force on the net charge Q 0 due to the uniform field 
E 0 . Hence 


F 

e 



q o E o 


(27) 
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in the Y direction, where Y = -R cos 0. Similarly, F h may also be calcula- 
ted by using a theorem (ref. 11) valid when Re << 1 and may be expressed in 
dimensional form along T direction as 


F 


h 


= -8^1 


1 im 1 
R->® R sin 


-e L 


+ - 


R sin 0 


(28) 


This completes the specification of the boundary value problem. It may be 
noted that the normal stress balance equation at the droplet surface has not 
been used in the formulation thus far. This is because the droplet is assumed 
to be spherical. Strictly speaking, it is the normal stress balance that 
determines the true shape of the droplet as it moves. We will return to this 
point later. 


SOLUTION 

(a) Case 1: Consider a case where q 0 2 /k2M1 << 1 as is the case for 
highly conducting droplets. The surface current due to charge convection may 
then be neglected and the problem is strictly electrostatic. Equation (25) 
reduces to 


8<j>, 

_£ 

8r 


0 at r = 1 


The solution to equations (13) and (14) is 


<!>, = -r cos 0 + 


H o 1 cos 0 


*,'o 7 * r 2 


^2 = Ci = constant 


= 


r ® 2 1 , , . 1 

- x + ^ — r + «( 1 - v ) - 
2 2 2 oo r 


sin 2 0 


= t (3v - 2)(r 4 - r 2 )sin 2 0 

C 4 00 


V = 
oo 


, .a 

h. 


1 + 


3p' 


2 Hi 


The dimensional terminal velocity is 


V = 


q o E o R o 




1 + 


P 2 


1 + 


3)i, 

2 Mi 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 


(35) 


which is the same as in reference 4. 

When + °°* there is no flow in the droplet and the results reduce to 
the motion of a solid sphere in a fluid. The terminal velocity of the solid 
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sphere may be derived by equating the Coulombic force on it to the Stokes 
drag. This yields 4irR 0 2q 0 E 0 = 6-irRQmV^, = 2q 0 E 0 R 0 /< 3pi ) , which checks with 

equation (35) in the limit p2 * 00 • From the solution, since <j >2 = and 
3^/30 = 0 at r = 1, there are no tangential Maxwell stresses and no surface 

tension gradients at the surface of the droplet. All the driving force for 
the flow comes from the normal Maxwell stresses, which is what contributes to 
F e . Since the normal stress condition at r = 1 has not been considered, the 
driving force is not explicitly seen. To see the effects of surface tension 
gradients, which is the primary goal of the present study, another case must 
be considered where charge convection must be retained. 

(b) Case 2: Here q 0 2 / k 2l i l not small and charge convection is 

retained. The resulting nonuniform charge distribution will induce surface 
tension gradients (cf, the molecular picture). However, for mathematical sim- 
plicity, the surface charge distribution is taken to be uniform in the charge 
conservation boundary condition, i.e., the change in the surface charge density 
is much smaller than the average surface charge density q 0 and 3q/30 is 
negligible. The exact condition, when this is true, is given following 
equation (40). For this case q * 1 and equation (25) becomes 


3<t>, 


3r ” k 2 p 1 sin 0 30 


Tx(vsin 0) at r = 1 


(36) 


The solution is 


♦ i - - r cos 0 + — r - + 

1 Vo r 


1 - 


Vi 


(3v - 2) 


COS 0 


(37) 


<j> 0 = - z (3v - 2)r cos 0 

T 2 k 2 p ] 


(38) 


v =- 

00 


3 + 


2(3q 0 3p 2 


3(3q Q 9p 2 

k 2 p 1 + 2p 1 


(39) 


y and <y 2 are the same as in equations (32) and (33), with v^ as given 
above. The surface charge density calculated from equations (21) and (37) to 
(39) is 


q = 



k 2 p 1 [ ill 

L ^2 

* k 2 H| + 2 M| 


COS 0 


(40) 
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Thus q « 1 is typically the case when c]E 0 , e2 E o << do- ma y be n °ted 
that in the limit q Q 2/(k2m) -*• 0, equation (39) becomes equation (34) and the 
solution for Case 1 is recovered. Also v® ■* 2/3 when p2 1 * °°. which is the 
result for a solid sphere. It is interesting to note that v„ also becomes 
2/3 when q 0 2 /(k 2 jii) -* 00 • This may physically be explained for the case with 
k 2 ■* 0 as follows. In equations (25) or (36), when k 2 -» 0, the conduction 
current at the droplet surface tends to zero, as we expect finite electric 
fields within the droplet. Hence, the divergence of the current due to charge 
convection, the charge convection current iteself, and hence the surface veloc- 
ity must be zero. This is precisely the case for a nonconducting solid sphere, 
where the surface velocity is zero on account of the no slip condition. In 
Case 2, driving forces due to normal and tangential Maxwell stresses at the 
surface and surface tension gradients are all present. As mentioned before, 
the effects of the surface tension gradients are contained in the tagged vari- 
able 13. 13 » 2 includes surface tension gradients, while (3 = 1 does not 

((3 = 1 may also represent a case where surface tension gradients are retained 
while Maxwell stresses are not; however, the Maxwell stresses must then be 
neglected everywhere in the formulation and would result in F e E 0 in the 
consideration of the net force on the droplet, a result that is not physically 
sensible). From equation (39), 


3v 


oo 


3|3 “ 



3q 0 2 3p 2 

+ + 2 h 


(41) 


Thus, for any |3, 3Voo/3(3 is always negative. Hence v® is smaller for |3 = 2 
than for (3=1, and so all other things held constant, inclusion of surface 
tension gradients results in a lower terminal velocity of migration of the 
droplet. The dimensional terminal velocity including surface tension gradi- 
ents is 



, V 3 ^2 

-Q 

O 

m 

O 

73 

O 

[ + ^1 + hj 


, 6q 0 2 9>1 2 


(42) 


(c) Deformations of the droplet: 

Having determined the terminal velocity and the velocity and potential 
fields, the perturbations to the shape of the droplet may now be determined, 
assuming that such deformations are small. As mentioned before, it is the nor- 
mal stress balance at the droplet surface that determines its true shape and 
will now be considered. The following analysis closely resembles that in 
references 13 and 14. The true shape of the droplet is represented by 
R(0) = R 0 + f(0). The normal stress condition may be written as 


N 


1 



at R = R + f 
o 


(43) 


where N is the normal stress due to both hydrodynamic and electrical effects 
and R a , R^ are the principal radii of curvature of the interface. When 
n = f/R 0 << 1, the following may be written (ref. 13) 
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/ e \ df 

^Re + t r e) de 


(44) 


N ‘ t RR + t RR 


R o + f 


1 

R 


a 


+ 





1 

R 2 sin 0 
o 


d_ 

de 


[»•« • I] 


(45) 


where trr = -P + 2jj 3U/8R, -c e RR = c/2(Er 2 _ E 0 2 ), P being the hydrodynamic 
pressure. The dimensionless pressure p = P/(jiVr/R 0 ) in the two fluids may be 
shown to be (ref. 1) 


P 1 = - ~2 cos 0, p 2 = K q - 20Aj r cos 0 (46) 

where A and A] are the coefficients of r sin 2 © and r 4 sin 2 © terms of the 
respective stream functions and K 0 is an unknown constant. The surface ten- 
sion a is a function of 0 and from equation (8), it may be written as 

o(0) = o Q - Y q [$(R q ,0) - (47) 

where q has been assumed to be independent of 0. y is again a tagged vari- 
able, where y = 1 retains the surface tension variations, while y = 0 regards 
the surface tension as uniform. 

Droplet shape for Case 1: Assuming n << 1 and using equations (44) to 
(47) and the solutions to (j> and y for Case 1, equation (43) may finally be 
written as 


l d_ 
sin © d0 


[sin 


0 

w de] 


+ 2 n = 2 - 


^2 V F 


K o " 


2e \°o 


9 C,E ° R ° rr- 2 n 

o COS 0 

2 cr 


(48) 


where the left hand side of equation (48) is essentially the right hand side 
of equation (43). Neither (3 nor y appear in equation (48) as 3<|>/8© = 0 
at r = 1 for Case 1. The boundary conditions on n are that the volume of 
the droplet is unchanged as it deforms and the origin of the coordinates is at 
the center of mass of the droplet. For n << 1, these become 




n sin 0 d© 


- °’ Jo n 


cos 0 sin 0 d0 = 0 


(49) 


From the analysis by Brignell (ref. 15), the solution is 

•? e 1 ? 
n (6) = i — ° (3 cos^ ©-I) 

O 


(50) 


Hence, n << 1 is valid when eiE 0 2 R 0 /o 0 << 1. The shape of the droplet is an 
ellipsoid of revolution, with an axis of rotation along the flow direction. 

To leading order, the ratio of the lengths of axes in the flow direction to 
that in the normal direction is 
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X = [1 + n<0°>]/n + nC90°>] ~ 1 + 


(51) 


9e, E^ 2 R^ 
1 o o 

8a. 


One sees, therefore, the droplet elongates in the flow direction. 


Droplet shape for Case 2: Similarly, using the solutions for Case 2 and 
assuming that q 0 E 0 R 0 /o 0 << 1 (i.e., A o/o 0 << 1) and n << 1, equation (43) may 
finally be written for the combinations (3 = 2, y = 1 or (i = 1 , y = 0 as 


1 d_ 
sin 0 d© 

where 


[stn 9 H * 2n - 2 - ^ K 


Ci e 2 r 


0 - [j„ + (j, - j„)cos 2 el 


°o 1X0 " 2a o rz T "" l " " 2 


(52) 


J T =9 + 12 7-*- 
1 k 2 Ml 


(2 - 3v ) 


4 - 


T 2 < 2 - 3v »> 2 

; 2 m, 


(53) 


J 2 - 


o / f\ . \ 2 


J k 


2 2 <2 " 3 V 

2 P 1 


(54) 


It is interesting to note that for both the combinations of (3 and y, no cos 
0 terms are present in the right hand side of equation (52). y does not 
appear anywhere else in equations (52) to (54) when Aa/a 0 << 1. The presence 
of a cos 0 term causes a logarithmic singularity in the solution for n at 
0 = 0,ir. Hence, both the formulations with and without surface tension gradi 
ents predict finite deformations of the droplet. The solution is 


n(0) 


1_ 

24 


g i E o 2r o 
1 0 o 


(Jl 


J 2 ) ( 3 cos 2 


0 - 1 ) 


(55) 


Once again, n << 1=^ E 0 2 R 0 /a 0 << 1 • However, q 0 E 0 R 0 /a 0 << 1 that has been 
assumed is a more restrictive condition as e i E 0 /q 0 << 1 for Case 2. It may 

be shown that in the limit q 0 2 /k 2 Pi -> 0, equation (55) becomes equation (50) 
and the droplet shape for Case 1 is recovered, thus serving as a check. Also., 


a eiE 2 R„ 
on _ 1 o o 

9(3 - 24 a Q 


( 3cos 2 0 - 1) 




3p 2 
+ 2^ 


(k 2M] ) 


123 - 10 + 


3d, 


^i 


3u, 

2 Hi 



(56) 
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and 3|n|/3(i > 0 for 3 > 5/6 - e2/(3ei>. Hence, |n| for |3 = 2 is greater 
than | n | for 3=1. Thus the inclusion of surface tension gradients increases 
the deformations of the droplet. Physically, this is to be expected as surface 
tension variations would tend to deform the droplet further away from its 
spherical shape. Under the conditions Aa/cr 0 << 1 and n << 1 , the increase 
in deformations comes from a greater imbalance in the radial Maxwell stresses 
at the surface of the droplet, originating from the inclusion of surface ten- 
sion variations. To leading order, the ratio of lengths of axes of the droplet 
in the flow direction to that in the normal direction is 


X 


e,E 2 R 

+ -4r 2(J i - V 


(57) 


DISCUSSION 

It has been assumed that the charge exists as a surface charge in the 
droplet. This is strictly true only when there is no flow and may be assumed 
to be approximately true when the convection of charge is small compared to 
charge diffusion in the bulk of the droplet. That is, charge conduction 
restores surface charge much faster than surface convection disrupts it. From 
the charge continuity equation in the bulk (ref. 10), this may be shown to be 
valid when e2VR/(k2Ro>» i-e*. PoEo^^m ^ • i- e -» Qo^^Pl) x e 2Eo/Po << 
Since ^o^Po must t> e small to justify the assumption of uniform surface 
charge distribution (eq. (40)), this imposes a limit on how big q 0 2 /(k 2 Pi) 
can be. 

Figure 2 shows the terminal velocity v (eq. (39)) for various values 
of ^/Pl and Po^^m) w ^th P * 2. As can be seen, v^ is of order one 
for all cases thus indicating that the reference scale chosen in appropriate. 
Also, v^ tends to 2/3, the value for a solid sphere, for large values of 

l^/pi or q 0 2 /(k- 2 Pi > . 

Figure 3 shows the shape of the droplet (eq.(55)) for v^Ml * 0-1* 
e l E 0 2 r o /c7 o b and various values of q 0 2 /(k 2 m) and zilt}. When 
q 0 2 /(k 2 ni) = 0, the shape is independent of ciU\, which is consistent with 
equation (50), the result for Case 1. The droplet can be elongated or con- 
tracted in the flow direction, the latter being the case for large values of 
tilzy and q 0 2 /(k 2 m>. The perturbations from a spherical shape must however 
be small for the theory to be valid. Many of the assumptions made are true 
for small values of E 0 ; hence, this theory is essentially a small field 
strength theory. 


SUMMARY 

The aim of this study is to examine the effects of surface tension varia- 
tions due to electrical effects on the migration of droplets in an insulating 
medium under microgravity condition. The important results are that surface 
tension gradients exist only when charge convection is included and tends to 
decrease the migration velocity and increase the extent of deformations of the 
droplet. The assumptions that the Reynolds number of the motion is small com- 
pared to one and that the charge redistribution is negligible essentially 
restricts the theory to small applied field strengths. Further analysis must 
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be performed to include the effects of electric double layers at the droplet 
surface, non-negl igible Reynolds numbers and the redistribution of charge on 
the migration velocity and shape of the droplet. 


REFERENCES 

1. Levich, V.G., Physicochemical Hydrodynamics , Prentice Hall, Englewood 
Cliffs, NJ, 1962. 

2. Levine, S. and O'Brien, R.N., "A Theory of Electrophoresis of Charged 
Mercury Drops in Aqueous Electrolyte Solution," Journal of Colloid and 
Interface Science , Vol . 43, No. 3, June 1973, pp. 616-629. 

3. O'Brien, R.W. and White, L.R., "Electrophoretic Mobility of a Spherical 
Colloidal Particle," Journal of the Chemical Society, Faraday Transactions 
II, Vol. 74, No. 9, 1978, pp. 1607-1626. 

4. Spertell, R.B. and Saville, D.A., "The Roles of Electrohydrodynamic 
Phenomena in the Motion of Drops and Bubbles," Proceedings of the 
International Colloquium on Drops and Bubbles , Vol. 1, D.J. Collins, M.S. 
Plesset, and M.M. Saffren, eds., Jet Propulsion Lab, California Institute 
of Technology, Pasadena, CA, 1974, pp. 106-121. 

5. Saville, D.A., "El ectroki neti c Effects with Small Particles," Annual 
Review of Fluid Mechanics , Vol. 9, M. Van Dyke and J.V. Wehausen, eds., 
Annual Reviews Inc., Palo Alto, CA, 1977, pp. 321-337. 

6. Hunter, R.J., Zeta Potential in Colloid Science - Principles and 
Appl i cations , Academic Press, New York, 1981. 

7. Callen, H.B., Thermodynamics , John Wiley and Sons Inc., New York, 1960. 

8. Hiemenz, P.C., Principles of Colloid and Surface Chemistry , Marcel Dekker 
Inc. , New York, 1977. 

9. Ono, S. and Kondo, S., "Molecular Theory of Surface Tension in Liquids," 
Handbuck der Physik , Vol. 10, S. Flugge, ed.. Springer Verlag, Berlin, 
1960, pp. 134-280. 

10. Melcher, J.R. and Taylor, G.I., "Electrohydrodyanmics: A Review of the 

Role of Interfacial Shear Stresses," Annual Review of Fluid Mechanics , 

Vol. 1, W.R. Sears and M. Van Dyke, eds., Annual Reviews Inc., Palo Alto, 
CA, 1969, pp. 111-146. 

11. Happel , J. and Brenner, H., Low Reynolds Number Hydrodynamics , Noordhoff 
Inti. Publishing, Leyden, Netherlands, 1973. 

12. Jackson, J.D., Classical Electrodynamics , John Wiley & Sons Inc., 

New York, 1962. 

13. Taylor, T.D. and Acrivos, A., "On the Deformation and Drag of a Falling 
Viscous Drop at Low Reynolds Number," Journal of Fluid Mechanics , Vol. 18, 
Mar. 1964, pp. 466-476. 


15 


14. Balasubramanlam, R. and Chai , A.T., "Thermocapillary Migration of 

Droplets: An Exact Solution for Small Marangoni Numbers," Journal of 

Colloid and Interface Science, Vol . 119, Oct. 1987, pp. 531-538. 

15. Brignell, A.S., "The Deformation of a Liquid Drop at Small Reynolds 
Number," Quarterly Journal of Mechanics and Applied Mathematics . Vol. 26, 
Pt. 1, Feb, 1973, pp. 99-107. 



FIGURE 1. - SKETCH OF A MIGRATING DROPLET. 
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